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The Lindblad approach to continuous quantum measurements is applied to a system composed of 
a two-level atom interacting with a stationary quantized electromagnetic field through a dispersive 
■ coupling fulfilling quantum nondemolition criteria. Two schemes of measurements are examined. 

' The first one consists in measuring the atomic electric dipole, which indirectly allows one to infer 

the photon distribution inside the cavity. The second one schematizes a measurement of photon mo- 
mentum, which permits to describe the atomic level distribution. Decoherence of the corresponding 
reduced density matrices is studied in detail for both cases, and its relationship to recent experiments 
is finally discussed. 
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Interaction between atoms and photons has been a fundamental issue since the early days of quantum mechanics, 
and it continues to be a central topic expecially in connection with controlled manipulation of small numbers of photons 
: and atoms at the quantum level In this framework, emphasis has been put recently on repeated measurements 
^\ ' on atoms and photons in a quantum regime requiring that the measurement process be explicitly taken into account. 
\ It is often desirable that the outcome of a measurement is not influenced by the previous ones, which has been 
found to be possible by a clever choice of the measured and measuring systems. Indeed, a class of measurements 
aimed at repeatedly monitoring the same observable has been introduced and shown to be compatible with the 
foundations of quantum mechanics, the so-called quantum nondemolition (QND) measurements Many QND 
' schemes have been proposed, and some of them have been implemented, to monitor various measurable quantities, 
^ , e.g. displacements of macroscopic oscillators below the standard quantum limit ^ , photon number in traveling jj] or 
O^' standing |^ electromagnetic fields, magnetic flux in superconducting interference devices vibrational energy of an 
^ \ atom confined in a Paul trap . 

From the theoretical point of view, the problem of incorporating quantum nondemolition measurements within 
the language of modern quantum measurement theory was not addressed specifically, to our knowledge, although 



^ . important steps were taken in |^,^. In this article, we analyze a model for QND measurements based on Lindblad 
■ ■ ' formalism, specializing it to measurement strategies aimed at describing continuous quantum nondemolition counting 
of photons confined in a cavity or atoms interrogated to be in a certain internal state. The system is described by 
a density matrix with an evolution equation introduced for generic open quantum systems pO| and already applied 
to other, demolitive measurement schemes such as the ones involving the quantum Zeno effect in hyperfine atomic 
spectroscopy [0, optogravitational cavities [Q, superconducting circuits ||T^, and trapped ions Q. In both cases 
the evolution of the coherence of the monitored system is studied, providing a simple picture of its decay during 
continuous QND measurements. Even starting from factorized states for atoms and photons, the measurement 
originates an entanglement that is ultimately responsible for the indirect decoherence of the observed system through 
the continuous collapse of the state of the probe. The implications of the model are discussed in connection to recent 
experiments implementing nondemolitive counting of photons in cavity QED jHj and atoms confined in electromagnetic 
traps 
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II. QUANTUM NONDEMOLITION COUPLING AND OPEN QUANTUM SYSTEMS 



A. General formalism 

In this Section we recall some useful concepts related to quantum nondemolition measurements, referring the reader 
to 1^ for a more detailed account. Although some measurements may occur through direct observation of the same 
observable under study, QND measurements are the most important example of indirect measurements, whereby the 
interaction of the system S with another system P is required. The dynamics of the monitored observable As of the 
system S is inferred through the modifications induced on the actually measured observable Ap of the probe. The 
total Hamiltonian operator is written as 

H^Hs + Hp + H^nt , (1) 
and the evolution of the overall system is described by the density matrix equation 

= . (2) 

The key idea behind a QND measurement is that the subdynamics of the observable Ag of the monitored system, 
albeit influencing the evolution of the probe, is not affected by this last. This nonrcciprocity is obtained if the 
interaction Hamiltonian Hint depends on As but docs not commute with Ap, i.e. 

[As,H^nt]^0, [Ap,H,nt]^0. (3) 

The subdynamics for the observables As and Ap are then ruled, in the Hcisenberg picture, by the equations: 

= [Hs, As] , = [Hp + H^nt,Ap] . (4) 

In addition, the observable As should be protected against the evolution of all the other system observables whose 
subdynamics is instead unpredictably affected by the measurement, in compliance with the Heisenberg principle. This 
requirement restricts the class of the observables that may be monitored in a QND way, a further sufficient condition 
being that the observable is a constant of the motion in the absence of interaction, i.e. [As, Hs] — 0. 

This description is still not a measurement model. In particular, the unitary evolution (Q) is in contrast with 
the general expectation that a measurement process should introduce irreversible signatures into the evolution of 
the system. Actually, the interaction whereby the system S is monitored by the probe P does not define alone a 
measurement process. As a next step the observable Ap has to be registered and, in order a measurement is defined 
univocally, this should correspond to a deterministic, classical, amount of information. It is the irreversibility implied 
by this act, with the pointer choosing only one of the many possible quantum alternatives, that cannot be accounted 
for within the closed dynamics discussed so far. At the beginning of quantum theory, this problem was solved by 
introducing a postulate, which in the von Neumann form states the instantaneous collapse of the wavefunction on 
the eigenstate corresponding to the observed eigenvalue [|6_ 
to various situations (including the cavity QED case as in 
a unified dynamics for quantum systems is built, without additional postulates, using the theory of open quantum 
systems. In this framework, the irreversible nature of the measurement is recovered by imagining the act of reading 
the outcome as due to an external environment which continuously interacts with the observed system. From a 
physical viewpoint, any process in which the interrogations on the observed system are made at a repetition rate 
larger than its intrinsic characteristic frequencies may be schematized as a continuous measurement. By restricting 
to so-called nonselective ensemble evolutions that can be represented in terms of a density operator, the most general 
dynamical law preserving (complete) positivity and normalization of the density operator was derived for a system in 
interaction with a Markovian environment in the form of a so-called Lindblad master equation [lO| ] , later successfully 
exploited in modeling quantum optics experiments [ p^ . Decoherence is introduced into this picture of open system 
evolutions as the dynamical quenching of the off-diagonal density matrix elements, a key property first used to explain 
the absence of superposition states in a measurement apparatus in [Q. Within this perspective, measured systems 
are nothing but a special class of open quantum systems, the environment being represented by the many-mode field 
of the macroscopic measuring apparatus and the measurement process being equivalent to repeated instantaneous 
effect- valued measurements or decoherentization kicks given to the density matrix We are thus led to 

the following master equation for the density of the coupled S + P system undergoing a measurement through Ap: 



This approach, although applied even in recent times 
H), has been overcome by a different description where 
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= ~l [ff (0, - ; Ap, \Ap,p{t)\ . (5) 
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The probe observable, Ap, plays the role of a Lindblad operator representing the influence of the external environment 
and the parameter k, with dimensions [k] = [i^^A^^], gives the couphng (generally time-dependent) of the probe to 
the measurement apparatus. By choosing a continuous function of time K{t) a continuous measurement process is 
obtained. 

B. Atom-photon Hamiltonian 

In the following, we will discuss quantum nondemolition counting schemes for both photons confined in a cavity 
and two-level atoms in a given eigenstate. By denoting with a, the standard annihilation and creation bosonic 
operators, the Hamiltonian of the free single-mode electromagnetic field is written as 

H photon = fiLoa^a , (6) 
while the two-level Hamiltonian can be expressed as 

Hatom — ^UJef(Tz (7) 

with respect to the basis of the energy eigenstates |e), |/), ojef — {Ee — Ef)/fi. It is convenient to shift the zero of 
energy in order to have Ef = 0, which is equivalent to choose the following shifted two-level Hamiltonian 

Hatom = fli^ef + ^z^ = TjCJe/CT+lT- = flLOeffle , (8) 

where (7+,(T_ and Hg — a^a^ denote the Pauli displacement operators and the projector over the excited state |e) 
respectively. Based on the previous considerations, we choose an interaction Hamiltonian which is linear in both the 
photon number operator a^a and the atomic projector Hg, 

H,nt = 2h-fa^alle , (9) 

the coefficient 27, a measurement angular frequency, quantifying the strength of the quantum nondemolition coupling. 
The occupation probability of level |e) for an atom in a generic state p("-tom) jg ^ Tr{p(''*°™)ne}. We note that, 
provided the Hamiltonian ^ is reinterpreted as single-particle operator for an ensemble of independent atoms and 
collective effects due to quantum statistics or interatomic forces are neglected, the probability Pg is related directly to 
the average number rie of atoms in level |e) through = PcTIt, nx being the total number of atoms. Once chosen as 
system observable Ag, both the photon number operator and the atomic occupation probability automatically satisfy 
the commutation relationship with the respective Hamiltonians (^), (H) since, in the absence of interactions, they are 
conserved. The total Hamiltonian to be used in is written explicitly as 

H = fiU)a)a + huJef^+<7- + '2h'ya) aa+a^ . (10) 

We will be working in the representation of the unperturbed (field -|- atom) eigenstates, expanding the density operator 
as 

00 

pii)= X! X! Pan,bm.{t)\an){hm\ , (11) 

a^b—e,f n,m— 

where 

Pan,bm{i) — (a7i|p(t) I^TTT-). It may be worth to look, for a moment, at the closed evolution of the density matrix 
elements. By projecting Eq. (||) (or Eq. (H) with k = 0), we find 

Pfnjm. = -iLu{n- m)pfnjm, 

Pen,e7n = -j(w + 2-f){n - m)pen,em, (12) 
Pfn,e7n = -i[w(n - m) - UJ^f - 2jm]pfn,em , 

whose solutions are simply rotations of the initial density matrix elements: 

Pfnjm{t) = exp [-iuj{n - ■m)t]pfnjm{0), 

Pen,eni{t) = CXp [-i{uJ + 2j){n ~ ■m)t]pen,eni{0), (13) 
Pfn,eni{t) = exp{-i[w(n - to) - UJef - 2jm]t}pfn^emi0) , 
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and Penjm{t) = P*fjn en(0- These equations imply obviously that both the average photon number and the average 
occupation probabilities of level |e) are time-independent, consistent with the QND nature of the interaction. On the 
other hand, the photon-atom interaction induces additional phase-shifts in the density matrix, which are proportional 
to the strength of the coupling 7 and depend upon the atomic and the photonic state. These phase-shifts contain 
the useful information about the system that needs to be extracted through the measurement on the probe. In 
the following two Sections, we will examine two complementary measurement procedures based on this optoatomic 
coupling. 



III. QND ELECTROMAGNETIC FIELD MEASUREMENTS VIA ATOMIC DIPOLE QUADRATURE 

A first class of QND measurements in optoatomic systems is obtained by monitoring the photon field using non- 
resonant atoms as the probe system. This corresponds to 



Hs — Hpho 



Hp — Hatom, -^s — a'' a. 



(14) 



Accordingly, the Hamiltonian (^0|) represents the fact that the photons in the cavity induce a state-dependent dy- 
namical Stark effect on the atom and a selective phase-shift of the atomic wavefunction. Due to the requirement 
of non-commutativity between the probe observable and the interaction Hamiltonian, we discard any operator pro- 
portional to the Pauli matrix tr^ as atomic probe operator Ap. On the other hand, in order to model a mea- 
surement which is sensitive to the dephasing accumulated between the components e and / of the atom inter- 
acting with the mode, phase-sensitive observables like the quadrature components of the atomic dipole operator 
ay — ((T+ — (7_)/2z, are natural candidates. Indeed, we observe that for a proper atomic superposition state, 
IV') = o|e) + ^1/); o, = |a|exp(z0a), b = |6| exp (i(/)b), |a|, |6| ^ 0, the average value (ay) = |a| |5| sin((/)a — </)&) is 
nonzero whenever a relative phase {(pa — 4>b) is present. By comparison, {a)x — \a\\b\ cos(0a — (j>h), maintaining finite 
values even if the relative phase vanishes. We therefore choose 



Ap^ 



0-4 



(7- 



2i 



(15) 



as the probe observable and by using Eq.(^ we get the following equations of motion: 

K 

Pfnjm = -iuj{n - m)pfnjr, 



^{Pfn.fm Pen, em); 



-i{uj + 2j){n - m)p 

en,e 



" {.Pen, em P/n,/m ) i 



(16) 



Pfn.em = -i[c^(n - m) - UJef " 2-fm]pfn,e 



^ (,Pfn,em ~1~ Pen, fro) 



Penjm = Pfrn en' uotc that the cvolutions for atomic diagonal and non-diagonal entries are decoupled but, at 
variance with the situations analyzed in pTj-p^, the measurement affects all components. 
By introducing the two families of pseudofrequencies u and w defined as 



u{n, m) — Y j^^n — m)^ 
Eqs. dig) are exactly solved giving 



16' 



w{n, m) — \ [ujef + "/{n + m)Y 



16' 



(17) 



PfnJ-ni{t) 
Pen, em 
Pfn,em (0 



-z(n — m) {ij+'y)t — Kt / 4 



_ g-i(n-m)([»j+7)t — Kt/4 



cos ut + i{n — m)j 
cos ut — i{n — to)7 



sin ut 
u 

sin ut 

u 



Pfnjni{0) 



K sin ut 



cos wt + i[LUef + {n + m)7 



pen, em (0) 

sin wt 



4 u 
K sin ut 
4 u 



Pen , em (0) 



P/n,/m(0) 



(18) 



w 



K sin wt 

Pfn,em['J) — ^ — Pen,fm['J) 



where the dependence of frequencies u, w on photon numbers n, m is understood. As expected, Eqs. ( |18| ) reduce to 
( p^ ) when no measurement is performed, k — 0, and for open-system coupling k small compared to 7, the dynamics of 
the coupled S + P system is just weakly perturbed with respect to the closed case. In the opposite regime where k 3> 7 
and a proper measurement on the signal mode is performed, pseudofrequencies u, w tend to become purely imaginary. 



4 



introducing overdamped oscillations and thereby decoherence. The time development of the electromagnetic field 
under the effect of the measurement can be inspected by evaluating the reduced density matrix: 



+ [cos^t + ( - z(n - m)7 + 9 ^] pif°")(0) ipH-^'Ho) 



K \ sin ut 



(19) 



where an initially uncorrelated state Pan,bm{0) = pl^i^°™\o)pnm^"'' (0) has been assumed. It is immediate to recognize 



(field) , 



that, for every k, field populations are unaffected by the measurement, i.e. pnn'^^'^\t) = pkn'^'"'' (0), in agreement with 
the QND nature of the coupling. However, the process of acquiring information on the field photon number cannot 
be realized without a back-action on the conjugate field variable, specifically an unavoidable degradation of the field 
phase distribution. A direct visualization of the phase evolution is provided by suitable phase-coherence indicators. 
We adopt here the formalism of the so-called Pegg-Barnctt phase distribution discussed at length in p3] : 



(field) I 



lim — 

s^oo 27r 



E 



(field) 
rnrn 



{t)e 



-i{n—rn)0 



(20) 



We have analyzed in detail the evolution arising when the field is initially in a coherent state jae*"^) with a mean 
photon number and the probe state is an equal superposition of levels e,/. As it will become clear in Section V, 
these choices will enable us a straightforward comparison with the results reported in jsj. We note, however, that 
equivalent physical insight would be gained from using a different standard quantum-optical distribution, the so-called 
Q- function Starting as a sharply peaked function centered sA 9 — (j), the Pegg-Barnett distribution retains its 
form when 7 = 0, the only changes reflecting the rotation of the coherent state in phase-space. When 7 7^ but k = 
the phase distribution is split into two components moving at different velocities with relative weights proportional 
to the coefficients of the atomic wavefunction. The evolution is strikingly different in the presence of measurement, 
since the phase distribution is progressively scrambled until it becomes completely flat when field coherences are 
asymptotically destroyed, pnm^'^\t) ^ 0,n 7^ m. The probability distribution ( pO| ) is shown in Fig. 1 for various 
instants of time and two different values of n corresponding to small and strong coupling with the environment. In 
the case of weak coupling (solid line), the splitting of the initial, well-defined phase distribution, into two components 
traveling with different velocities is still recognizable. 



IV. QND ATOMIC MEASUREMENTS VIA PHOTON MOMENTUM 

A second class of QND measurements is obtained if, according to the notations of Section IIA, we choose 

Hs — Hatomi Hp = Hphoton, = ^e- (21) 

This corresponds to monitoring the atomic level via an electromagnetic field and is realized ordinarily by means of 
an absorption process, i.e. by sending photons resonantly tuned at an energy-level gap. This technique is manifestly 
demolitive for the interrogated atoms. An alternative approach consists in the detection of the phase-shift originated 
by the atom through a non-resonant interaction with a light beam, the atomic sample acting as a dispersive medium 
with a complex refraction index. A nice application of this technique has been recently reported in ||l5|| , where 
repeated, nondestructive optical imaging of an atomic cloud confined in a magnetic trap has been demonstrated. The 
modification of the refraction index, proportional to the atomic optical density, induces a phase-shift of the probe field 
which is manifested in turn as a change in the amplitude of its quadrature component. The dispersive phase-shifts 
may be also detected through amplitude modulation by means of classical, standard techniques like dark-ground or 
phase-contrast imaging [ p5[ . In our model, it is natural to choose the following phase-sensitive photon operator, 

(22) 



2i ' 

as a probe observable. This choice, in addition to fulfilling the QND criteria established above, corresponds to a 
measurement configuration which is the dual of the one analyzed in Section IV. By exploiting the master equation 
(^ again, and by evaluating the new Lindblad commutator, equations of motion of the following form are derived: 
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Pan,bm — Pan,bm 



(k=0) _ K 
4 



{n + m+ l)pan,b'm + \/ n{m + l)pan-l,bm+l + \J m{n + l)pan+l,bm-l 
nnipan-l,brn-l ^ V (" + + ^)Pan+l,bm+l - ^ ^i^ - l)Pan-2,6m (23) 



+ a/ m{m - l)pan,bm-2 + \/ {u + l){n + 2)pan+2,bm + \/(m + l)(m + 2)pan,bm+2 

where a,b — e, f and n,m > 0, and p]^~^^ indicates for brevity the appropriate unmeasured contribution, shown 

exphcitly in Eqs. ([T^). Since, as before, the probe observable is diagonal in the system variables, no transitions are 
induced by the measurement process in the system. In this case, the e and / components evolve independently, but 
a quite complicated structure of the couplings in the Fock space of the photons is present in general. Eqs. (|2^) have 
been integrated numerically for a field coherent state and an equal-weight atomic superposition state. The reduced 
density matrix that is appropriate to study the evolution of the signal mode during the measurement is the atomic 
density obtained by tracing the total density matrix over the photonic degrees of freedom. 



pT'^\t) = Y.P'^-.bn(t)- (24) 
n=0 



By analogy with the previous case, we expect that although no perturbation is introduced on the variable which 
is QND-monitored, i.e. pit^°"^\t) = pie*°™''(0), atomic coherence is eroded due to back-action. In Fig. 2 the time 
dependence of the coherence indicator |pe/(i)| is plotted for different values of the coupling parameter k. The evolution 
of atomic coherence in the closed case k = can be evaluated explicitly for a coherent probe state since, from Eqs. 

1 °° 1 

Pef(t) = - e*""^* ^ e^*''"*pi{i*''''''' (0) = ga^(cos27t-l) ^i{ui^ft+a^ sin 27*) _ (^25) 

This shows that the quantity |pe/(i)|^ oscillates with angular frequency 27 between a maximum value equal to 1/4 
and a strictly positive minimum equal to e^^" /4, leading to an example of nontrivial reversible dynamics of the 
modulus of the atomic coherence. In the presence of measurement, atomic coherence damps exponentially to zero 
with a rate proportional to the parameter k, still preserving the oscillatory behaviour visible in Fig. 2. 



V. EXPERIMENTAL COMPARISON 



The model developed here, describing decoherence induced by a continuous QND measurement, is characterized 
completely by the parameters 7 and k. In this Section we discuss to what extent these parameters can be related to 
a realistic experimental scheme. The example we take is the scheme proposed and implemented by a group at the 
Ecole Normale Superieure where atoms detect in a nondemolitive way the photons stored in a high-Q cavity. In 
their scheme, the atom is schematized as a system with three Rydberg levels e, /, and i, with levels i and / of the 
same parity, opposite to the one of level e. The frequency lo of the cavity mode is detuned from the e — > i transition 
by an amount 6 = to — \ujie\- The detuning is large enough to neglect photon absorption, but small in comparison to 
the angular frequency . The presence of photons in the cavity results in a phase-shift of the e-state relative to the 
/-state. We recall that the dynamical frequency shift induced on an atom in level e and atomic dipole d and located 
at point r in a cavity containing N photons is 

Ae(r,iV) = ^| 

where E{r) is the electric field at r and the last approximation holds if absorption processes are made negligible. The 
phase-shift is proportional to the vacuum Rabi angular frequency fi(r) = E{r)d/h. Noticing that N is the eigenvalue 
of the d^d operator, we can think of the spatially averaged phase-shift as due to the effective Hamiltonian 

Hm - h '' ; " d^dU,. (27) 




1 + V N 



E{ryd' 



N = 



nir) 



-N 



(26) 
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The parameter 7 of our model is therefore identified as 7 = (J7^(r))/25. Concerning the parameter k, the discussion 
is more elaborated, since in the realistic situation analyzed in |^ the measuring meter is actually a beam of two-level 
atoms crossing the cavity and the effect of the measurement is taken into account using the von Neumann collapse for 
each atomic interaction. Distinction is made between two different configurations of the probe system, corresponding 
to either atoms with a known initial velocity that are registered by the field ionization counters, or atoms which are 
not in a monokinetic velocity state and are not read. In the language of quantum measurement theory, this difference 
can be restated in terms of selective and nonselective measurements nonselective dynamics being obtained by 
averaging over all the possible states of the probe, in this case over the velocity distribution measurements and the 
two possible outcomes of the internal level. To make explicit comparison with Q , it has been shown there that the 
selective evolution of the field density matrix at the {k + 1)*'' atomic detection event is given by: 



iT\a,v) = ^^S^^^^^pitUa,-), (28) 

Z^n I'^ay'lT', V)\ 



where ba{n, v) denotes the amplitude of the component a {a — e, f) corresponding to a generic atomic velocity v and 
a photon number n. As first step towards nonselective measurements, we need to evaluate the weighted average over 
the possible final outcome of the probe, leading to: 

pI^^'Hv)- E &a(n,.)6:(m,.)p«(.) = {sin^^+cos^^e-("-™)^}pW(.), (29) 

a=ej 

where the parameter e measures the accumulated phase-shift per photon and vq is the atomic velocity corresponding 
to a 7r/2 pulse in the Ramsey zone used for interferometric detection of the dephasing. The next step is a second 
integration over the atomic velocity distribution P(v): 

pit'' = Jdv Piv) E ba{n,v)b:im,v)pl^iiv). (30) 

a=ej 



In Figs. 3 and 4 the Pegg-Barnett phase distribution is plotted for a monokinetic atomic beam, (Eq. (|2S|)) and a 
thermal atomic beam, (Eq. (^o|)), respectively. While in both situations the phase distribution is broadened and 
tends towards a flat one, only in the second case the transient does resemble the behaviour found for a nonselective 
measurement by using the effective Lindblad approach, Cfr. Fig. 1. In particular, the velocity of decoherence is 
proportional to the temperature of the atomic beam, i.e. the beam variance. This qualitatively agrees with the 
result established in p| , concerning the direct proportionality between the measurement coupling constant and the 
temperature of the bath in which the meter is embedded [ pT| . Unlike the parameter 7, it is not possible to infer a 
simple relationship relating k to the various parameters of the realistic configuration. This shows advantages and 
disanvantages of the effective Lindblad approach: it is often impossible to relate it completely with the experimental 
set-up; however, a general description of any measurement process is obtained without a detailed knowledge of the 
actual experimental procedure. For this reason, the formalism may be easily adapted to the description of other 
relevant schemes, like the single photon-atom coupling in a high-finesse Fabry-Perot cavity ||2^, the QND counting 
of atoms in an optogravitational cavity based upon use of evanescent fields p5| , p6| , and the nondestructive imaging 
of a Bose-Einstein atomic condensate . For the latter system, the dephasing of the atomic coherence induced by 
the measurement process should generate damping of the oscillatory behaviour of either a two-species condensate 
exhibiting Rabi-like oscillations or two single-component condensates spatially separated and undergoing coherent, 
Josephson-like oscillations, a situation similar to the one involving a superconducting circuit already analyzed in [[l3[. 



VI. CONCLUSIONS 



A model for the description of the decay of coherence in continuous quantum nondemolition measurements involving 
photons and atoms has been developed and applied to photon and atom counting. Even if initially uncorrelatcd, atoms 
and photons become entangled during the measurement and, although conserving their number, they are subjected 
to decoherence due to the back-action via the probe. The dynamics of the system has been analytically solved in the 
case of photon counting via atomic detection. The decay of atomic coherence has been evidenced when the atomic 
occupation probability is measured through the monitoring of the light beam. Contact has been established with actual 
experiments in which decoherence induced by the measurement process is or could be observable. In particular, the 
nondestructive monitoring of Bose-Einstein condensates of atomic dilute gas through dispersive imaging, demonstrated 
in will deserve a particular attention, since the influence of the measurement on the dynamics of the phase of 
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the wavefunction is a crucial issue in the study of macroscopic quantum coherence. More in general, this model gives 
constraints on the minimum rate at which coherence is progressively destroyed during a continuous measurement 
process even if a quantum nondemolition monitoring is adopted, provided that all other possible sources of decoherence 
have been quenched by proper technological improvements. As a consequence, similar considerations should be also 



relevant to investigate decoherence dynamics within QED-based quantum computation proposals 28 1 or quantum 
control strategies involving QND-mediated feedback Finally, we point out that in the model described here 

emphasis has been put on considering the whole coupled (system+probe) object as an open quantum system in 
interaction with a measuring environment, the probe degrees of freedom being traced out on the total density matrix 
and not dynamically eliminated from the beginning. This makes the description formally different from other existing 
approaches MM . A quantitative discussion of the equivalence between the two strategies will be addressed elsewhere 
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FIG. 1. Time evolution of the Pegg-Barnett phase distribution for a coherent state of the electromagnetic field and various 
values of the measurement coupling constant k — 10^^ (solid line) and k = 10 (dashed line) in the continuous case. The first 
20 energy eigenstates have been used as a truncated basis for the photon field, leading to a numerical accuracy of 0.1%. The 
snapshots from above to below are taken at times differing by one measurement period defined as Tm = 277/27. 

FIG. 2. Square modulus of the atomic reduced density matrix versus time (in units of the measurement period Tm) 

for various values of k. The oscillations are present also in the closed case {k = 0), although not visible in the scale, and 
corresponds to the behaviour of Eq. (25). 

FIG. 3. Time evolution of the Pegg-Barnett phase distribution for a Von Neumann measurement on a coherent state of the 
electromagnetic field (with initial Pegg-Barnett phase distribution as in Fig. 1), corresponding to consecutive interrogations 
with a monokinetic beam at two different velocities, v/vo = 0.3 (solid) and v/vo = 0.6 (dashed). The phase is progressively 
scrambled although its spreading starts from well-defined regions. 

FIG. 4. Time evolution of the Pegg-Barnett phase distribution for a Von Neumann measurement on a coherent state of 
the electromagnetic field (with initial Pegg-Barnett phase distribution as in Fig. 1) with thermal atoms for two different 
temperatures (in arbitrary units) T = 10~^ (solid) and T = 10 (dashed). The scrambling of the phase is faster than in the 
monokinetic case, affecting all the phases, and is proportional to the temperature. 
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